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In pressure sensors, four piezoresistors connected in a Wheatstone bridge are 
often provided with a voltage varying between 5 and 10 V. Unfortunately, 
this voltage is a source of drifts created by electric heating. This study 
focuses on the internal heating and piezoresistive effect of piezoresistors 
represented by the variation of their resistivity. To do this, we use the finite 


difference method (FDM) to solve the heat transfer equation, taking into 
account the conduction in Cartesian coordinates for the variable regime. We 
Keywords: examine how the temperature affects the piezoresistivity in these sensors 
when the potential is applied. In this case, the variation of the temperature 
has been calculated as a function of applied voltage, as well as for the 


Finite difference method 


Mobility sue operating time of the sensor. Furthermore, the evolution of resistivity over 
Piezoresistivity time was determined for several geometric properties of the membrane using 
Sensors the mobility model. This was established for different doping levels. 
Voltage Additionally, the change in resistivity due to the application of voltage was 
evaluated. It was observed that resistivity is greatly affected by the 
temperature rise produced by the applied voltage when the device is actuated 
for a prolonged time. Consequently, this results in drifting in the output 
response of the sensor. 
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1. INTRODUCTION 

Many domains that require the use of piezoresistive sensors with excellent reliability such as 
industry and biomedical due to their good linearity of output voltage, excellent compatibility with 
microelectronic systems and some other advantages [1]-[4]. However, they often suffer from the thermal 
drift. The study of electric heater effect in these sensors types is important to optimize the output 
characteristics drift [5]—[8]. As a result, several previous research papers have been carried out in this domain 
in order to lessen this effect, when the bridge is powered by an electric tension. Research by 
Ansari and Cho [7] was to derives both temperature independent and dependent models for temperature rise 
created by bias voltage in piezoresistive microcantilevers. Their approach is based on analytical model in 
steady state for describing the temperature diffusion in microcantilever by self heating [7]. Research by 
Aryafar et al. [9] introduce a novel technique for temperature compensation of piezoresistive pressure 
sensors. In their technique, additional silicon resistors with negative temperature coefficient of resistivity are 
used for compensation purpose. More recently, Tian et al. [10] have developed an analytical model of 
temperature characteristics in pressure sensor to investigate the influence factors of temperature zero drift, 
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including temperature dependence of resistivity. In their paper, they have studied the influence of several 
piezoresistor size combinations on the zero point output of pressure sensors using the finite element method 
(FEM) technique to effectively compensate for temperature zero drift [10]. To highlight the effect of 
electrical heating in piezoresistive pressure sensor, Pramanik et al. [11] have established an analytical model 
for self heating in circular-shaped diaphragm sensors. Nevertheless, the sensitivity of the gauge to the 
constraints of the square form diaphragms sensors is greater than in circular shape when the diaphragm 
thickness is greater than 2 um [4], [12], [13]. Similarly, at the same geometric properties of membrane 
(length or diameter) the circular-shape diaphragms has greater rigidity compared to those of square shape 
[14]-[16]. As a consequence, square form membranes are more sensitive than circular ones [4]. For this 
reason and some others uses such as fields that require excellent sensitivity to pressure, the diaphragms of 
sensors are usually made using square or rectangular shape [16]-[18]. This prompted us to present a 
numerical model of Joule heating in piezoresistive pressure sensors, with a square shaped membrane 
[19], [20]. We have analyzed the effect of supply voltage, geometric properties of membrane and the 
piezoresistor length, this allows us to improve the performance of this sensors. Recently, in the next paper 
[21] we have studied the impact of thermal drift on sensitivity of pressure sensors. The obtained results allow 
us to optimize the pressure sensitivity of sensors by optimization the geometric properties. Considering the 
importance of piezoresistivity effect in the performance of piezoresistive sensor a study was carried out by 
Ansari and Gangadhara [22]. In their paper, a finite element analysis is used to find the design parameters of 
the piezoresistor for application in atomic force microscope and some other applications. Accordingly, 
Kerrour et al. [23], have explored the change of resistivity due to temperature for several doping 
concentrations based on the mobility model of Arora. Their approach, has made it possible to optimize the 
response of piezoresistive pressure sensors of thermal drifts when used in high temperature areas. So as, to 
study the reliability of piezoresistive pressure sensors, we have coupled the numerical model of self heating 
with that of the mobility of Arora to optimize the influence of bias voltage on the piezoresistivity effect, 
described by the change of the resistivity of the piezoresistor, when applying mechanical stresses. Based on 
the numerical resolution of 2D heat transfer equation using FDM, we examine the geometrical, electrical and 
physical parameters influencing the resistivity to optimize the thermal drifts in the response of this sensors. 
For that, the evolution of temperature due to the self heating was found as a function of operating time of 
sensor and of applied voltage. Moreover, we have calculated the variation of resistivity in time for several 
geometric properties of the membrane and for different doping concentrations, as well as, for bias voltage. 
Afterwards, we embarked on the study of the influence of temperature on sensitivity. 


2. THEORY AND GOVERNING EQUATION 

Piezoresistive pressure sensor is composed of four gauges placed on the membrane in order to get a 
high piezoresistivity; thus, to get excellent sensitivity. Thereby, two gauges are placed parallel and the two 
others perpendicular to the membrane edges as shown in Figure 1. These gauges are connected in a 
Wheatstone bridge. 


Pin Metal base Narrow beam Petal membrane 


Wheatstone bridge ——. 


Piezoresistors —— 


Center boss Pyre 7740 glass 


Figure 1. Structure of piezoresistive pressure sensor [14] 


The piezoresistivity effect is strongly influenced by the temperature; particularly, that generated by 
the self-heating, which pushed us to concentrate on this source of heating. This involves the resolution of 
heat transfer equation in cartesian coordinates for the transient regime. This equation is governed by [7]: 
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where k is thermal conductivity, Q is heating flux, p is the mass density, t is the time and C is the specific 
heat. Otherwise, the rate of energy production by electric heater is given in (2) [7]: 


Apzr Vo 
= 2 
Q Lpzr-d-Sm Pe ( ) 


where Vo is the bias voltage, d is the diaphragm thickness, Sm = a? is the diaphragm area, Lpzr is the length 
of the piezoresistance, Apzr the cross-sectional area and p, is the electrical resistivity expressed in (3) [23]: 


_ 1 
Up (T)Na 


Pe (3) 


where q is the elementary charge, Up is the mobility and N4 is the doping concentrations. The starting 
condition in the whole structure as in (4): 


Tia, y,t=0)=T (4) 


The boundary conditions including the adiabatic heat condition and preserve the heat continuity at the bords 
can be given as in (5): 


ar _ ðr 
— = — = 0 
eee, (5) 


After discretization in two dimensions of the heat transfer equation in transient state using FDM, we find a 
system of linear equations which is solved by Thomas method using the MATLAB software. Within the 
framework of study of the reliability of these sensors, we are interested to the thermal behaviour of resistivity 
during their operation time. So, we use the formula of mobility established by Arora [24], [25]: 


p(T) = Um — 1s (6) 


1+(74) 


where Umn is minimal of mobility, uo, the difference between the maximum and minimum value of 
mobility, Nen the concentration reference and 0 is an exponential factor are given by [25]: 


estes eae? (7) 
Hon = 1250 ee (8) 
Nan = 1.26 x 1027 je)" (9) 
6 = 0.88 a" (10) 


By substituting the mobility expression (6) in the resistivity (3), we obtain the final formula of the resistivity 
as a function of temperature and doping level. 


pe = —> (11) 
on Hon sJ 


NA 
1+(74) 
Nen 


3. RESULTS AND DISCUSSION 
As we have noted in the previous section, the thermal drifts induced by electric heater in 
piezoresistive pressure sensors affect highly the response of such sensors. This work aims to investigate the 
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impact of the thermal drifts in this sensors when, bias voltage between 5 and 10 V is applied. It also, put 
emphasis in the study of the geometric influence parameters and doping level on the piezoresistivity effect. 


3.1. Influence of bias voltage and operating time on temperature rises 

In this section, the study of effect of bias voltage and of operating time of sensors has been 
performed. Firstly, we have represented in the curve Figure 2, the variations of temperature provoked by the 
bias voltage established by the numerical model at t=100 min and that of the analytics [7]. The comparison 
analysis shows a good correlation between the analytical and numerical model (FDM) for the low voltage. 
Indeed, beyond 8 V there are some differences, which may be due to the mesh, whom it is a source of error in 
the numerical model. From this figure, we have found that, the higher the bias voltage is, the greater the 
temperature rise will be, which is in accordance with the theory. 


32 F F F T 


= Numerical model 
== Analytical model [7] 


2 E É lad 
65 7 8 
Voltage (V) 


Figure 2. Effect of bias voltage on the rise of temperature 


Secondly, as reported in Figure 3, it can be observed that the change in temperature induced by self 
heating has an exponential form. It is noted that, from t=100 min we are in permanent mode and that the 
evolution of the temperature stabilizes becoming independent of time. Additionally, the shape of this curve is 
analogous to the response of a first order circuit. The latter, is characterized by a time constant of the order of 
3T from which the steady state is achieved. It is described as in (12): 


POSE. ses) (12) 


From (12), we can make a correspondence between the temperature T and the potential Ec. Therefore, we can 
get a relationship between T and the operating time of the sensor. 


T(t) =Tn(1—e7) +T, (13) 


where Tn and T are the constant and the thermal time constant. 

Thirdly, the drifts of resistivity when the piezoresistive pressure sensor is in a high temperature 
environment was investigated by [26], [27]. In our case, we have plotted in Figure 4 the variation of 
resistivity versus temperature rise created by the electric heater, when the bridge is powered by an electric 
tension. We can notice a significant increase in resistivity according to the temperature induced by the 
self-heating. Thereby, drifts in characteristics of the sensor. 
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VO=5 V, a=1000 um, d=20 um, Lpzr=40 um, Apzr=6 um2 


F 


E É sal 
100 150 200 
Time (min) 


Figure 3. Effect of operating time in the rise of temperature 
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Figure 4. Variation of resistivity versus temperature rise generated by electric heater 


3.2. Influence of time and geometrical parameter of membrane on the resistivity pe 

The knowledge of resistivity variation when the stress is applied is used to quantify the evolution in 
the output voltage, produced by these constraints. By coupling the temperature variation in time from the 
used model in (13) with the Arora mobility (6 and 11), we obtain the change in resistivity pẹ in operating 
time of sensors for the p-type silicon (Na=4.10!° cm-3). The evolution of the resistivity in time and in 
geometric parameters of membrane are illustrated in Figures 5-7. From the Figure 5, we can observe that the 
Pe is an increasing function with time. This increase is most certainly due to self heating. It is further 
observed in Figures 6 and 7 that greater the thickness and wide of the membrane are, lesser is the resistivity 
variation. However, to increase the pressure sensitivity, thickness of the silicon membrane must be reduced. 

According to these Figures 6 and 7, the two parameters have an tremendous effect on the resistivity 
and, hence, on the piezoresistivity. Thus, when these geometric parameters are great, this leads to reduce the 
self heating. In contrast, it is important to remember that this solution can have consequences on the 
mechanical and technological behaviour of the sensors, specifically on its pressure sensitivity and reliability 
on the one hand and on its size on the other hand. 


Bulletin of Electr Eng & Inf, Vol. 12, No. 5, October 2023: 2784-2793 


Bulletin of Electr Eng & Inf ISSN: 2302-9285 O 2789 
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Figure 5. Resistivity variation with operating time 
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Figure 6. Resistivity variation vs membrane thicknesses 
Na=4e16 cm-3, VO=5 V, d=20 um, Lpzr=40 um, Apzr=6 um2 
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Figure 7. Variation of resistivity in time for different membrane wide a 
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3.3. Influence of time and of doping concentrations on the resistivity 

In the same way, we have varied the operating time and doping level with the fixing of the other 
parameters. We found that the resistivity decreases by increasing the doping level as shown in Figure 8. 
Afterwards, the p, becomes independent of the operating time for high doping concentrations. So, the change 
in resistivity can be reduced by increasing the doping concentration. 
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Figure 8. Resistivity variation vs time for different doping concentrations 


3.4. Influence of applied voltage on the resistivity 

The effect of the bias voltage in the resistivity is plotted in Figure 9. It is noticed that the pẹ observes 
an increase when the applied voltage is increased. This is apparently due to the phenomenon of internal 
heating in the piezoresistor. As we have pointed out earlier, the applied voltage involves a rise in temperature 
subsequently it leads to an increase in resistivity. To summarize, the increase in resistivity due to the 
application of a voltage of a bridge appeared low, but it leads to drifts in the sensor response. 
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Figure 9. Resistivity variation vs voltage 


3.5. Effect of temperature rise generated by electric heater on sensitivity 
The output characteristics of the piezoresistive pressure sensors in terms of sensitivity is of great 
importance. The sensitivity of these sensors is an essential characteristic to optimize their output response. In 
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the center of the square-shaped diaphragm of thickness d and length a, their sensitivity is expressed 
in (14) [28]: 


S(T) = 0.1479 a T44(No, 300K)P(N,T) (14) 


where 744 is the piezoresistive coefficient established by [29]. 

As mentioned earlier in Figure 3, for a bias voltage of 5 V and for 100 min of the operating time of 
the device, the temperature rise caused by electric heater can reach 45 °C. This encouraged us to study the 
influence of temperature on sensitivity when the temperatures rise from room temperature to 45 °C. 
Accordingly, it can be established by coupling the change of temperature as a function of time from the used 
model in (13) with (14). It has also been observed in Figure 10 that the sensitivity decreases almost linearly 
with temperature. Subseqeuently, it causes a drift in the output characteristics of the sensor upon application 
of a voltage. So, the results revealed that the self-heating effect on the sensitivity is significantly reduced for 
a short operating time and for a low applied voltage. 
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Figure 10. Variation of sensitivity vs temperature rise induced by self-heating 


4. CONCLUSION 

In summary, this research has confirmed that it is possible to model, accurately the thermal drift 
provoked by bias voltage in piezoresistive pressure sensors, by coupling the finite difference model for a self 
heating whit these of the Arora mobility. The model developed gives an opportunity to study the temperature 
rise created by supply voltage in such sensors. The results proved that low bias voltage should be applied to 
lessen electric heater on the piezoresistivity effect, which is characterized by change of resistivity. They 
furthermore show that the evolution of temperature is proportional to the time and takes a steady state value 
from 100 min. So, the self heating effect is reduced substantially for a short operating time of sensor. 

In addition, we have attempted to study the geometric influence parameters on the resistivity to 
optimize the performance of the device. The results showed that the resistivity decreases as the geometric 
parameters of the membrane increase. Therefore, self heating can be reduced, when these parameters are 
great. Nevertheless, these parameters are themselves limited by: the dimensions of the sensor, the pressure 
sensitivity and reliability. It can be also reduced as a high doping level as possible it is used. 

To sum up, this paper allows us to evaluate the furthermore. We have also endeavored to study the 
influence of the temperature rise induced by the Joule heating on the sensitivity of the sensor. The results 
confirmed that the sensitivity is a decreasing function of temperature rise. So, to reduce thermal sensitivity, a 
low supply voltage must be applied. Reliability of the device it equally permits to improve the output 
characteristics of the sensor according to the application for which it is dedicated. 
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